Astronomy & Astrophysics manuscript no. annihil'DM'CMB 


© ESO 2013 


January 14, 2013 





o 

CD 

a^ 

(N 



O 

u 

o 



(N 

cn 
o 



- '—I 

X 
J3 



WMAP7 and future CMB constraints on annihilating dark matter: 

implications for GeV-scale WIMPs 

Gert Jens Chluba^, Andi Hektor^, Martti Raidal"^ 

^ Tartu Observatory, Toravere 61602, Estonia 
e-mail: gertSaai . ee 

^ Max-Planck-Iristitut fiir Astropliysik, Karl-Scliwarzsdiild-Str. 1, D-85741 Garching, Germany 
^ Canadian Institute for Theoretical Astrophysics, 60 St. George Street, Toronto, ON M5S 3H8, Canada 
e-mail: jchlubaOcita.utoronto . ca 

National Institute of Chemical Physics and Biophysics, Tallinn 10143, Estonia 
e-mail: andi .hektorOcern. ch, majrtti .raidalOcern. ch 



Received / Accepted 

ABSTRACT 

Context. We calculate constraints from current and future cosmic microwave background (CMB) measurements on 
annihilating dark matter (DM) with masses below the electroweak scale: moM = 5 — 100 GeV. In particular, we assume 
the S-wave annihilation mode to be dominant, and focus our attention on the lower end of this mass range, as DM 
particles with masses moM ~ 10 GeV have recently been claimed to be consistent with the CoGeNT and DAM A/LIBRA 
results, while also providing viable DM candidates to explain the measurements of Fermi and WMAP haze. 
Aims. We study the model (in)dependence of the CMB power spectra on particle physics DM models, large-scale 
structure formation and cosmological uncertainties. We attempt to find a simple and practical recipe for estimating 
current and future CMB bounds on a broad class of DM annihilation models. 

Methods. We use a model-independent description for DM annihilation into a wide set of Standard Model particles 
simulated by PYTHIA Monte Carlo. Our Markov chain Monte Carlo calculations used for finding model constraints 
involve realistic CMB likelihoods and assume a standard 6-parameter ACDM background cosmological model, which 
is extended by two additional DM annihilation parameters: muM and {(Jav) /moM- 

Results. We show that in the studied DM mass range the CMB signal of DM annihilations is independent of the details of 
large-scale structure formation, distribution, and profile of DM halos and other cosmological uncertainties. All particle 
physics models of DM annihilation can be described with only one parameter, the fraction of energy carried away by 
neutrinos in DM annihilation. As the main result we provide a simple and rather generic fitting formula for calculating 
CMB constraints on the annihilation cross section of light WIMPs. We show that thermal relic DM in the CoGeNT, 
DAM A/LIBRA favored mass range is in a serious conflict with present CMB data for the annihilation channels with 
few neutrinos, and will definitely be tested by the Planck mission for all possible DM annihilation channels. Also, our 
findings strongly disfavor the claim that thermal relic DM annihilations with m-oM ~ 10 GeV and {(jav) ~ 9 x 10~^^ 
cm'^s^^ could be a cause of Fermi and WMAP haze. 

Key words. Cosmology: theory - dark matter - diffuse radiation - cosmic microwave background - Elementary particles 

1. Introduction particle- and astrophysics communities to try to find other 

than purely gravitational manifestations of DM. So far, 

According to our current understanding, the energy-matter we have good knowledge of DM only through its gravita- 

content of the Universe is dominated by dark components: tional effects, starting from the scale of galaxies and galaxy 

dark energy (DE) 73%) and dark matter (DM) (- 22%), dusters, up to the cosmologically largest observable scale s 

with ordinary baryonic matter accounting for only ~ 5% (jJungman et al.lll996t iBertone et all l2005t lEinastol [2009t ) . 

of th e total (~ critical) density (e.g., see [Komatsu et_al] However, as there is already an impressive list of ongo- 

|201lD. In contrast to the poorly understood DE - the sub- ing and upcoming direct DM detection experiments along 

stance causing the Universe to expand in an accelerated with various indirect means of detection (see Pens ,2010l 

fashion - we have physically well- motivated models for the for overview), the hopes are very high that in the near- 

DM. Among those, the most promising scenario states that est future the mystery of DM might at last be solved, 

the DM of the Universe consists of thermal relic density Indeed, the first signals from DM particles could poten- 

of stable weakly interacting massive particles (WIMPs). It tially have aheady been detected: the (ex pected) annual 



ui nuduic wccuviy iiiuciaLuiiig, iiiaoBivc jjdi uilicb yvviivii b^. i.u tiaiiy uave aireaoy been aetected: tne (exp 
is quite miraculous that having particles with masses and modu lation signal from DAM A /LI BRA (iBernabei et al 



annihilation cross sections set by the electroweak scale au- |2010D . signals from the CDMS (| Ahmed et all 1201 ll 



tomatically provide the right DM density after freeze-out and CoGeNT (jAalseth et all 1201 ih nuclear recoil exper 

dJungman et al.||1996|; [Bertone et alJ|2005D. iments, anomalies of the cos mic ray positr o ns as re- 

The WIMP hypothesis, along with its potentially ob- vealed by PAMELA satellite (|Adriani et all l2009t) . or 

servable phenomenology, has initiated strong effort in the positrons-|-electrons as obtained by the Fermi satellite 
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(jAbdo et all 120091) and HESS atmospheric Cherenkov 
telescope (lAha ronian et all 120091 ) . While the cosmic ray 
positron an omaly can possibly be expla ined by T eV-scale 
DM (Bcrast rdmet al.ll2008l:lBarger et a l. 2009; Circl li et al] 
[2009k .Arkani-Hamed et alJ 120091: iFox fc Poppitz, 12009^ 
the signal from CoGeNT calls for liRht WIMPs withi n 
the mass ranee 7-12 GeV (lAa lseth et all l201lt) . 
Consistent analyses of combined data from CoGeNT and 
DAMA/LI BRA determine t he light DM mass to be 6 - 
8 GeV (|Hooper et all l2010D. Although this mass range 
is probed by the C DMS (|Ahmed et al][20TTI). XENONIO 
(Angle ct al. 2009^, and XENO NIOO jAprile^eLaU [20Tot 

experiments, inter- 



KENO.NIOO Collaboration et al. 2011 



pretation of those results ( Ko pp et al. 



^oTot iSchwetllmol) 



requires an ability to reliably reconstr uct nuclear recoils 
at very low energy ([Savage et al.ll201lh . as well as precise 
knowledge of DM distribution and velocity in the local halo. 
Therefore the CoGeNT and DAMA/LIBRA hints of light 
DM cannot be ruled out unambiguously. 

In addition, there is an independent positive claim of 
th e existence o f 0(10 ) GeV mass DM. A recent study 
bv lDobler et al.l ([201 1[ ) also suggests that annihilating DM 
with similarly low masses (ttidm = 1 — 20 GeV) may give 
a good match to the observed Fermi and W MAP haze 
([Dobler et al.l 120101: iDobler k Finkbeineil [2008') . However, 
all those claims depend strongly on the knowledge of the 
profile of the DM halo of our Galaxy and precise knowledge 
of local DM density and halo substructure. 

Thus the several interesting claims of the existence of 
0(10) GeV mass DM call for model-independent tests of 
the light DM scenario. Since lower DM particle masses 
imply higher number densities (riDM c>c 51dm^^/"^dm), and 
as the energy input from annihilations scales as oc n'^-j^, 
one might expect strong constraints on annihilation cross 
section, which might possibly reach below the standard 
thermal production v alue of {oav) ~ 3 x 10" 



-26 



([Jungman et al.|_ 19961). The co nstraii its from gamrna-ra; 

Bell & Jacques' '2009^ 'Bertone et a ll 1200' 
2009; Cirelh k Panci 2009; Meade et al 



measurements 
Bergstrom et al. 
20101: ICirelli et al.l 



_ 2009; 
20101 iPapucci Strumiai l2010t 



Hiitsi et all 120101: [Baxter et al.l l2010t lArina k TvtgatI 



20111: IVincent et al.H2010t IZavala et al.fl2011l) alon g with 



CMB b ounds (Galh et al.' '2009; 'Slatver et al.' '200; 
ICirelh et al. 2009a; Hiitsi ct al. 200i); Kanzaki et al. 201 
indicate that this might indeed be the case. 

In this paper we investigate how well DM annihilation 
cross sections for WIMPs with masses moM — 5 — 100 
GeV can be constrained with CMB measurements, in par- 
ticular focusing on the lower end of this mass range. The 
main advantage of CMB over other indirect probes, like 
gamma ray and cosmic ray measurements, is that it is prac- 
tically insensitive to the complications caused by the non- 
linear evolution of the cosmic density field. As our analy- 
sis shows, for any realistic structure formation scenario the 
CMB bounds on annihilating DM arise solely around the 
redshifts of z ~ 1000, while the contribution from lower 
redshift cosmic structures is completely negligible. 

CMB constraints on annihilating DM have 

been obtained in seve ral e arlier studies : e 

Padmanabhan &: Finkbeineil (I2005D : iMapelli et all ([20& 



Zhang et al.' (2006"); 
2009): Cirelli et aL 



'Gain et al.' 
(2009a) : 



2009"): ISlatver et alJ 



^ Hiitsi et al.l (|200E 

iKanzaki et al.l ( 2010() . Even though most of these analyses 



the energy depositioiFI , more recent studies ([Slatver et al.l 
120091 : iHiitsi et "al]l2009l ) followed the energy transport prob- 
lem including various energy-loss mechanis ms in a more 
realist ic way. Compared to the analysis of ISlatver et al.l 
([2009f ). which partially relies on t h e pre viously derived 
'on the spot' results of I G alii et al.l (|2009D . in this paper 
we perform a more elaborate treatment for the energy 
deposition joined to the analysis of the CMB data via 
Markov chain Monte Carlo calculations that incorporate 
the most recent WMAP likelihood code. Also, we make 
an attempt to unify the results from various annihilation 
channels and provide a simple and rather generic fitting 
formula for calculating CMB constraints on annihilation 
cross section (aAv) for a broad range of annihilating DM 
models. 

Our paper is organized as follows. In Section 2 we give 
a brief description of the energy input from DM annihila- 
tion and provide a simple treatment for its propagation. 
The effect on CMB temperature and polarization fluctu- 
ations is investigated in Section 3. Section 4 presents our 
main results about current and future CMB constraints. 
Our summary is given in Section 5. 



2. Input signals and their propagation 

As in ICirelli et ah! (|2009bl l201ll) . we treat our input sig- 
nals from DM ann ihilation in an as m odel independent a 
way as possible. In ICirelh et"all ([201 iD the two -particle an- 
nihilation channels to all Standard Model (SM) particles 
were considered: leptons, quarks, photons, gluons, weak- 
interaction gauge bosons, Higgs boson, and neutrinos. In 
addition, annihilations to four leptons via an intermediate 
new boson V were considered. Such a treatment can be 
considered as model independent since realistic models can 
always be decomposed into these basic channels where the 
particular branching ratios between the channels are given 
by the underlying theoretical particle physics model. 

Since in our work we focus on DM particle masses below 
100 GeV, out of all of the above channels the following re- 
main: DM DM SM SM, where SM= {e, /x, r, q, c, 6, 7, g}; 
plus 4-lepton channels via V. Here q denotes the light 
quarks u, d, and s. Because the masses of interest in this 
work are mostly below the masses of the electroweak gauge 
bosons, and Z, those channels are left out. For the same 
reason we also do not need to distinguish between left- and 
right-handed particles. Also, the neutrino channels in this 
case provide only trivial output; i.e. 100% of the energy 
is carried away by neutrinos, which escape freely at the 
redshifts of interest, so are not treated any further. Even 
though the channels 7 and g are included in our model- 
independent approach, these are strongly suppressed for 
realistic models since DM should not carry color or interact 
electromagnetically. 

For all channels, the spectra of the emerging stable par- 
ticles, e, p, 7, after treatment of several decays, parton 
showers, and ha dronization were calcu lated using PYTHIA 
Monte CarldEI (I Siostrand et al.l l2008[) . All input spectra 
are downloadable from |http: //www.marcocirelli .net/| 
,PPPC4DMID .html. For more details, and in particular for 



have assumed a simple 'on the spot' approximation for 



^ I.e., energy input from DM annihilations is assumed to get 
instantaneously absorbed by the cosmic medium, with some ef- 
ficiency factor /. 

^ |http : / /home . thep . lu . se/ ~torb j orn/Pythia . html 
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Fig. 1. The redshift z where the optical depth for photons 
reaches unity (i.e., t(z, z ) = 1) for several 'observer's red- 
shifts': z = 0, 10, 500, 1000. Here the energy plotted is the 
photon energy at redshift z. The light gray region corre- 
sponds to the DM mass interval considered in this paper. 



a discussion on the level of possi ble uncertainties, we refer 
the reader to lCirelh etlll (|201lD . 

Among the stable output particles, neutrinos propa- 
gate freely at the redshifts of interest, while immedi- 
ately interact with the ubiquitous CMB photons and up- 
scatter those to the gamma-ray energy range via the in- 
verse Compton (IC) mechanism. The total output energy 
in hadrons {p and d) is typically quite negligible. Only in 
gluon and quark channels does it reach up to 15%, and 
that only for the highest DM particle masses considered in 
this paper. We therefore did not model this component in 
detail. However, in Sect. |4] we give some arguments that at 
the redshifts of interest, i.e. z ^ 1000, probably the ma- 
jority of energy released in p and d channels is directly 
converted into heat by the cosmic medium. Thus, we only 
need a detailed treatment for the photons: (A) the ener- 
getic ones originating directly from the annihilation event 
(prompt photons), (B) the softer IC photons created by 
upscattering CMB photons. 

The processes and the corresponding cross sections 
relevant to the propagation of photons through the cos - 
mic medium were taken from'Zdziarski &: SvenssonI (jl989ll . 
Starting from the lowest of energies these include (i) pho- 
toionization, (ii) Compton losses (on both bound and free 
electrons), (iii) pair production on matter, (iv) photon- 
photon scattering, and (v) pair production on ambient pho- 
ton fields. 

In Fig. [1] we show the redshifts z where the optical 
depth t{z,z ) = 1 for various redshifts of the observer: 
z = 0, 10, 500, 1000. As can be seen, for intermediate photon 
energies (depending on the redshift) there is a well-known 
X-ray/gamma-ray energy window where the photons can 
pro pagate freely over cosinologic ally large distances (e.g., 
see IChen fc Kamionkowski|[2004D . For the cosmological ra- 
diation transfer it is crucial that this 'transparency window' 
is properly modeled. Once the photon gets outside of this 
window (we take it to happen after the first interaction) we 
assume that the following cascade will be locally absorbed 



in a very short time. Moreover, the fractions (1 — /ion)/3 
and (1 -I- 2/ion)/3 of the total absorbed energy, are as- 
sumed to be going for ionization and heating, respectively. 
Here /ion is the fraction of ionized hydrogen atoms, and 
a similar expression for helium can be used. Excitations of 
atoms are neglected. This approxima tion, motivated by the 
workof lShuilfc van Steenberd (1198 5*1 , has been widely used 
in seve ral subsequent pape rs, e.g. Chen & Kamionkowskl 
(2004); Padmanabhan & Finkbeiner (2005); Mapclh ct aT| 
|2006i); .Zhang et al.. (,2006) : ,Nataraian fc Schwarz, (2009( 1 
However, it is clear that these simple expressions only pro- 
vide a rough estimate for t he corre c t ener gy deposition ef- 
ficiencies. As mentioned bv lChlubal (|2010[ ). the precise red- 
shift dependence of the efficiency factors and the fraction 
of energy that goes into excitations of hydrogen and he- 
lium atoms need more careful consideration of the radiative 
transfer processes, including secondary low-energy photons 
and their feedback. These extra photons have the potential 
of further delaying recombination, hen ce affect the last scat - 
tering surface and CMB anisotropics (jPeebles et al.|[2000l ). 
We leave a more detailed investigation of these ambiguities 
to a future paper; however, later on in Section |4] we briefiy 
comment on how much th e omission of excitation s using 
the rough prescription of IChen fc Kamionkowskil ()2004[) 
changes our final results. 

In Fig. [1] we also see that at high redshifts, as the 
density of the environment becomes much higher, the X- 
ray/gamma-ray transparency window starts to close. Thus 
at sufficiently high z, we would expect all the produced an- 
nihilation energy (excluding the energy stored in neutrinos, 
since these can freely leak out) to be absorbed locally. In 
the following we call the ratio of the locally produced to 
the locally absorbed energy the /-parameter. At high red- 
shifts (but well after the neutrino decoupling) we expect the 
/-parameter to asymptote to the value given by (1 — f^), 
where /^ is the fraction of energy carried away by neutri- 
nos. An example for the /-parameters in the case of /i an- 
nihilation channel are shown in the upper lefthand panel of 
Fig. [21 Since ^ 60% of the energy is carried away by neutri- 
nos in the /i-channel the expected asymptotic high-redshift 
/-parameter should be ~ 0.4, which is indeed the case. 
With this in mind, we see that robust model-independent 
results from the CMB analyses can be obtained for the DM 
masses below ^ 100 GeV. This is the reason we concen- 
trate on light WIMPS in this work. For heavier WIMPs, 
the computation depends on more complicated details for 
energy absorption. 



3. Effect on CMB 

Because we are able to calculate /-parameters, we can go 
on to calculate the effect on cosmological recombination. 
To this end we modify the cosm ological recombination code 
RECFAST (Seager et al."l999') along the lines presented in 
Pad manabhan fc F inkbeiner (2005). For more details see 



also IHiitsi et al 



()2009l) . Although more advanced cosmo- 

l ogical recombination codes have recently been releasee^ 
(iChluba fc Thoma£l[20Tot lAli-Haimoud fc Hiratall20TTh . at 
this point we do not include any of the recently discuss 

^ CoSMOReC (www.Chluba.de/CosmoRec) already provides a 
simple module that accounts for the effect of DM annihilation 
or more general energy injection, but this module is still being 
validated. 
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Fig. 2. Lefthand column from top to bottom: (i) /-parameters for the /j,-channel assuming ttidm = 10 GeV (dark solid 
lines) and muM = 100 GeV (light solid lines). In both cases, the lowest curve out of the triple of lines corresponds to 
only the smooth background, the middle one includes halos with a lower mass cutoff of 10~^Mf:^, while th e top one has 
\0~^Mq- The dotted lines represent high-redshift fits (valid for z > 170) as given bv ISlatver etall ()2009[ ). The dashed 
line shows the CMB visibility function. Here zg{z)d\nz gives the probability that the CMB photon last scattered in 
the logarithmic redshift interval dlnz centered on redshift z. (ii) Fraction of free electrons as a function of redshift for 
niuM = 10 GeV and 100 GeV, along with two values for (aAv) / {{<TAv)std'mY)M[GeV]): 1 and 10. Here {aAv)std — 3x 10^^^ 
cm'^s"^ is the standard thermal cross section. For the meaning of each line, see the description given in the legend. The 
lowest dotted line represents the standard ACDM case with no additional energy input from the annihilating DM. (iii) 
Matter temperature as a function of redshift. The models shown are exactly the same as in the panel above. The long 
dashed line compares the temperature of CMB. (iv) The CMB visibility function g{z) times redshift z for the same 
models as shown in the above two panels. Righthand column from top to bottom: (i) Angular power spectra of 
CMB temperature fluctuations for the same models as already given above, (ii) Temperature and E-mode polarization 
cross-spectra, (iii) E-mode polarization spectra. In all of the righthand panels the points with errorbars show the 7-year 
measurements by the WMAP space mission. 
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corrections to the cosmological recombination process aris- 
ing from de tailed radiative trans fer and atomic physics 
fsee lFendt e' t al. 2009; Rubiiio-Ma rtm et al.ll20Tol and ref- 
erences therein for overview). 

The main effect of DM annihilations is a delay in recom- 
bination at 2: ~ 1100 and an increase in the low-redshift 
freeze-out tail. This changes the position and width of 
the last scattering surface and thus affects the CMB tem- 
perature and polarization anisotropies (see, for example, 
Chen fc Kamionkowskil[2004l : iPadmanabhan fc Finkbeined 
20051) . Since the input annihilation power scales as oc p|)ji^, 
it is clear that structure formation leads to a signifi- 
cant boost over the average density squared p\)m^ i-S- 
{p\,m{^)) — B{'^)PDMi^)- The details of how w e calculated 
the str ucture boost factors B{z) can be found in lHiitsi et al.l 
(j2009t) . The onset of structure formation at z ~ 100 is 
clearly visible in the upper lefthand panel of Fig. [21 Here 
the three dark (light) solid curves, representing the results 
for the n annihilation channel, correspond to TOdm = 10 
GeV (100 GeV). From the bottom to top, the trio of lines 
in each case represent: (i) no structure formation, i.e. only 
a smooth background, (ii) including structures with the 
lower halo mass cutoff of 1O~^M0, (iii) structures with the 
cutoff of 1O~^M0. The halo mass funct ion was assumed 
to ha ve an analytic Sheth-Tormen form feheth fc Tormeiil 
19991) and the halo mass-concentration relation followed 
Maccio et al.l (|2008[ ) d escription. For the full details of this 
calculations, again see lHiitsi et al.l ()2009l ). The dotted lines 
represent the high-redshift fits (valid fo r z > 170) for the 
/-parameters as given by 'Slatver et al.l (|2009l) . After ad- 
justing slightly the high-redshift normalizations we see that 
the shape of their /-parameters at high redshifts agrees re- 
markably well with our results, even though our treatment 
is somewhat more simplified. The sharply peaked dashed 
curve at z ^ 1000 in Fig. |2] shows th e Thomson visibil- 
ity function (jSunvaev fc Zeldovichlll97dil . More specifically, 
the quantity V ^ zg{z)d\n.z gives the probability that the 
CMB photon last scattered in the logarithmic redshift inter- 
val dlnz centered on redshift z. For the CMB calculations 
the values of the /-parameter matter only for the redshift 
range over which the visibility fu nction is large, and in this 
range our results agree with the ISlatver et al.l (|2009[ ) cal- 
culations to better than 5% accuracy. We also tested the 
results for the r-channel and found very good agreement. 
Thus, our somewhat mo re simplistic calculat ion of the /- 
parameters compared to ISlatver et al.l (|2009[ ) seems to be 
justified. 

The other lefthand panels of Fig. [5] (from top to bot- 
tom) show the fraction of free electrons, matter temper- 
ature, and the Thomson visibility function as a function 
of redshift for two DM particle masses, 10 GeV and 100 
GeV, along with annihilation cross sections {10, 100} times, 
and {100,1000} times the standard thermal cross section 
((T^i;)std = 3 X 10^^^ cm'^s"-^. See the legend for the line 
definitions. The dotted lines show the standard behavior of 
ionization fraction, matter temperature and visibility func- 
tion without any additional energy input from DM anni- 
hilation. In the panel for the matter temperature we also 
show the line for the temperature of the CMB. 

Because the annihilation power scales in proportion 
to ((TAw)nQ]y[mDM cx {<JAv)/m-aM, we would expect the 
case with ttidm = 10 GeV, and {oav) — 10((TAu)std 
(100(cryiz;)std) give comparable results to the case ttidm = 
100 GeV and (ctaw) = 100(f7Aw)std (lOOO(crAw)std), as long 



as the /-parameter does not depend strongly on todm- This 
is indeed approximately so, as can be seen from Fig. [2] Even 
though the lines for the ionization fraction, matter temper- 
ature, and CMB visibility function are easily separated at 
lower redshifts, this simple scaling with {aAv)/iTiuM seems 
to hold very well near the peak of the visibility function, and 
so one would not expect to be able to clearly distinguish the 
models with the same value of {(TAv)/m-DM via CMB mea- 
surements alon^l. Also, the additional annihilation boost 
from the structure formation is not expected to influence 
the CMB measu rements. Only in th e very extreme case 
where instead of lMacci"6 "eTall (l2008l) mass-concentration 
relation we use a simple power-law extrapolation down to 
the very low halo masses one is able to cause an addi- 
tional significantly high low-redshift peak in the visibility 
function. However, for any realistic mass-concentration re- 
lation, along with annihilation cross sections {(tav) that 
do not violate the CMB data, the contribution from the 
structure formation to the CMB signal is completely negli- 
gible. Also, the contribution to the low-redshift ionization 
fraction is very mild, and thus the annihilating DM models 
which are compatible with CMB measurements could only 
play a marginal role in reionizing the low-z Universe. These 
results supp ort our similar findings previously reported in 
iHiitsi et al.l (|2009f ). 

In the righthand panels of Fig. [2l we show the 
temperature-temperature (TT), temperature-E-mode po- 
larization (TE), and E-mode-E-mode (EE) angular power 
spectra. The points with error bars show the WMAP seven- 
year measurements (jLarson et al.ll201l[) . The dotted lines in 
all of the panels show theoretical predictions for the con- 
cordance ACDM model (Komatsu et al. 2011). As we might 
already expect, remembering the visibility function behav- 
ior as shown in the lowest lefthand panel, the models with 
the same values for {(Tav) /muM are hardly distinguishable. 
This is because the redshift dependence of the /-parameter 
over the width of the visibility function is small, and so 
only the average value around z ^ 1100 really matter s, in 
agreement with the statements of ISlatver et al.l ()2009l ). 

As a final note in this section we point out that, if one 
increases (aAv) to high enough values and adds the struc- 
ture formation boost, so that the low-redshift ionization 
fraction starts to get significantly approach one, the original 
RECFAST code, along with its DVERK solver, is not prop- 
erly able to deal with the underlying set of stiff ODEs and 
simply breaks down. To be able to numerically treat these 
cases, one can use the significantly improved RecfastH — h 
codtQ ([Chluba fc Thomas 2010). along with its stiff ODE 
solver. However, in our calculations the original RECFAST 
works fine, since the typically allowed values of {<jav) / ttidm 
compatible with the CMB data are low enough. 



^ To achieve better distinction between various models, one 
might try to use the differences in the behavior of the low-z 
matter temperature and its observable consequence on the 21 
cm transition measurements of the neutral hydrogen. The effect 
of annihilating DM on t he 21 cm signal has been discussed e.g. 
in JvIa poUi ct al. (200 6): lFurla.netto et aP (|2006l ): IValdes et all 
(|2007l ): lNataraian fc Schwarzl l|2009l )" 

^ RecfastH — |- is part of more advanced recombination code 
CosmoRec and can be downloaded from www . Chluba . de/ 
CosmoReCj 
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Fig. 3. The WMAP7 parameter constraints for the /i annihilation channel. Along with two annihilation parameters 
log]^Q(mDM) and (f7yiu)/TODM, we also show the constraints on the other 6-parameter ACDM background model: flBh^, 
^Buh'^, Hq, t, ng, log]^Q(As). From the inside out the colored 2D areas show the 1-sigma and 2-sigma regions after 
marginalization over the other parameters, and the dashed lines show the same regions for the basic 6-parameter ACDM 
without annihilating DM. The topmost panels in each column plot the marginalized ID probability distributions for all 
of the parameters. 



4. CMB constraints 

Now we have all the ingredients available to calculate cur- 
rent and future CMB constraints on annihilating DM. We 
notice that the CMB constraints for the decaying DM are 
not competitive (e.g. when compared to the bounds avail- 
able from the low-redshift gamma-ray measurements) due 
to (X pdm scaling when compared to the oc p%]\.j scaling for 
the annihilating DM, and thus are not being discussed in 
this work. 

For the currently existing CMB data, we used the lat- 
est power spect ra, i.e. seve n-year spectra, from the WMAP 
space mission (jLarson et al. 2011). In reality the current 
CMB bounds can be somewhat tightened if one includes 
results from other smaller scale CMB measurements, in par- 
ticular the measurem ents of the temperatur e power spec- 
trum from ACBAR (jReichardt et al.l l2009l ) and temper- 
ature and polarization spectra from QUaD (jPrvke et all 
I2009D In addition data from ACT ()Das et al.ll201lD and 
SPT (jLueker et al.ll2010l ) could be added to the analysis. 



providing additional leverage on small scalesQ. However, to 
keep our analysis clearer we decided to restrict ourselves 
to the WMAP7 data only. Regarding future CMB results 
we made predictions for the bounds available from the cur- 
rently ongoing Planck mission, and also for the idealized 
noise- free experiment up to th e multipole of ^max = 2000. 
As in IZaldarriaga et aIT(|2008[) . for Planck we assume 80% 
sky coverage, beam size ^fwhm = 7 arcmin, and noisebias 

iV/T ^ IQ-i ^EE ^ 3 5 X IQ-i ^^2 foj. rprp 

EE, respectively. For the ideal noise-free experiment, we as- 
sume a full sky coverage and the underlying uncertainties 
of the CMB fluctuations are solely due to a finite number 
of available fluctuation modes on the sky, i.e. due to cosmic 
variance, and so we call this type of idealized experiment 
'cosmic variance limited' (CVL) in the following. First steps 
towards CVL CMB measurements in both temperature and 
polarization down to small scales will become available from 



Since part of the effect of DM annihilation is degenerate with 
the effect of changing the spectral index of scalar perturbations, 
ns, small-scale CMB measurements help break this degeneracy. 
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a combinati ons of SPTPOlFI (iMcMahon et all I2OO90 and 
ACTPoiI!! (jNiemack et al.ll2010t) along with Planck data, 
so that the CVL case discussed here provides a good guide- 
line to what could become possible in the near future. 

We generated the synthetic data for the Planck and 
CVL experiments as described in http : //cosmocof f ee . 

0/0. Along with WMAP7 likelihood codcE], they will be 
used in our Markov chain Monte Carlo (MCMC) parame- 
ter estimation. For the MC MC engine we used the publicly 
available CosmoMC tooliFlf Lewis fc Bridld I2002D . where 
the standard recombination modules were modified in order 
to allow additional energy input from annihilating DM- Our 
baseline concordance ACDM model (jKomatsu et al.ir2011[ ) 
is parameterized by six free parameters: flBh'^, ^Duh'^, Hq, 
T, ns, log]^o(^s), where fie and fluM are density parameters 
for baryons and DM, Hq — 100 • h km/s/Mpc is the Hubble 
constant, r is the reionization optical depth, and As and 
ns are the amplitude and spectral index for adiabatic scalar 
perturbations, respectively. This minimal six-parameter set 
is extended by two additional parameters describing the an- 
nihilating DM: its mass moM, and thermally averaged an- 
nihilation cross section (cr^u). Since MCMC will give better 
results if model parameters are as uncorrelated as possible, 
we used the parameter pair {todm, {o'Av)/mY)M} in place 
of {mDM, {o-Av)}- 

As an example, we show the WMAP7 constraints on 
all of the eight parameters, assuming /i annihilation chan- 
nel on Fig. O With dashed lines we show parameter con- 
straints for the six-parameter ACDM model without anni- 
hilating DM. The colored 2D areas in all of the panels dis- 
play 1-sigma and 2-sigma regions marginalized over all of 
the other parameters. For all of the parameters the topmost 
panels in each column plot the marginalized ID probability 
distributions. It is reassuring to see that after introducing 
two additional parameters, all the previous six parameters, 
albeit with small shifts, are still as precisely determined. 
The strongest shifts compared to the baseline six-parameter 
model are quite understandably seen for parameters Ag and 
ns, which both lead to the rise of the CMB power spectra 
if increased, and thus counterbalance the damping of the 
spectrum caused by the extra scattering off the additional 
free electrons created by the energy input from the DM 
annihilation. 

Here it is important to note that these small shifts 
towards higher values of As and ns are directed in the 
same direction as recently discussed recombination cor- 
rection s from previously neglected standard physics pro - 
cesses (|Rubiiio-Martfn et al.li20Tot iShaw fc Chlubal [201ll) . 
Although biases in Ag and ng arising from recombination 
physi cs are not significant for WMAP7 (Shaw & Chlub^ 
1201 If ), parameters derived from Planck data will be affected 
by several standard deviations if recombination corrections 
are neglected. This indicates that precise constraints on 
models with annihilating DM should probably account for 
these recombination corrections simultaneously. 



http : //pole .uchicago . edu/ 
* http : //www. physics . princeton. edu/ act/ 
^ See also http://lpsc.in2p3.fr/perotto/ in case CMB 
lensing is needed. 

|http : // lambda. gsfc .nasa. gov/produ ct/map/ current/ 
likelihood_i iif o . cf m 

[http : //^sinologist . info/ cosmomc/ 
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Z 

Fig. 4. Upper panel: /-parameters for all of the channels 
considered in this work and for several DM particle masses 
in the interval m-DM = 5 — 100 GeV. For all of the cases, 
the asymptotic high-redshift /-parameter values have been 
renormalized to be equal to one. Lower panel: The same 
as above, after dividing with a 'typical shape', (/max(-z) + 
/min(-z))/2, of the /-parameter curve. 



In Fig.[3]we only present the constraints for one particu- 
lar channel. However, we would like to provide easily usable 
recipes to calculate bounds for all of the channels described 
in Section [2] As seen before, even though the /-parameters 
vary somewhat if one changes the mass from 10 GeV to 100 
GeV, the CMB spectra as shown in Fig. [5] were fairly insen- 
sitive to this change as long as the ratio {aAv)/mY)M was 
kept the same. In the upper panel of Fig. [j] we plot the /- 
parameters for all of the channels and for several masses in 
the range todm = 5 — 100 GeV. Here we have renormalized 
the high-z values to be equal to one. In the redshift range 
where the visibility function peaks (see the lowest lefthand 
panel of Fig. [5]), the variation among different DM models 
is not too large. To see this better, in the lower panel of 
Fig. |3]we have divided the above lines by a 'typical shape' 
of the /-parameter given by (/max(2) + /min(^))/2, where 
/max(2) and /min(^) are the extremal values of f{z) for a 
particular z. We see that, at most relevant redshifts, the 
variation among all the different models is around 15%. 
This indeed suggests that it might be possible to provide a 
unified approximate scheme for calculating the bounds for 
all of the annihilation channels. 

To bracket all the possibilities, we calculated the 
bounds assuming the /-parameter to have a typical shape 
(/max(-z) + /min(-z))/2 aloug with two cxtrcmal values 
fininiz) and fmax{z). Of coursc, instead of two annihila- 
tion parameters we then only have to deal with one ad- 
ditional parameter: ((TAw)/mDM- We took the asymptotic 
values for the /-parameter to have values of 0.25, 0.50, 0.75, 
and 1.00. This asymptotic value, as explained above, is sim- 
ply (1 — /^), where fi, is the fraction of the total energy 
carried away by neutrinos. It turned out that the resulting 
1-sigma and 2-sigma constraints on (1 — f^){aAv) /todm for 
the above four values of (1 — /^) are the same with (1 — 2)% 
accuracy. Our final results (averaged over these small vari- 
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Table 1. l-sigma and 2-signia values of r in Eq. ((T)) and 
their uncertainties for WMAP, Planck, and CVL CMB ex- 
periments. 





r 


l-sigma 


2-sigma 


WMAP 

Planck 

CVL 


+ U-Ui>i 
0.073_o 

+ 0.0055 
0.0160_o 0022 
„ +0.0020 
O.OOri_o 0013 


+ U.UBB 
0191-0 031 

+ 0.0096 
0.0326_o 0055 

+ 0.0031 
0.0137_o.o032 



ations due to changing (1 — f^)) can thus be given in the 
form 



, ^ iaAv) [3 X 10-26 cm^s-i] 
(1 ~ U) FTT-TTl < 



(1) 



where the values of r for l-sigma and 2-sigma bounds are 
given in Table [TJ One can see that the approximate bounds 
we provide are typically accurate at about the 20 — 30% 
level. We tested that, indeed, the ranges of r given in 
Table [T] fully cover the values of r for all of the channels con- 
sidered in this work. Two examples for the case of l-sigma 
upper bounds are shown in Fig. [S] Here the upper panel 
corresponds to the /i-channel {f^ ~ 0.61) and the lower one 
to the e-channel (/^ ~ 0.02). We chose the above two chan- 
nels because these are the two extreme cases among all 
of the channels treated in this paper. Consequently, they 
bracket the expected results for any realistic annihilating 
DM, given as a superposition over the basis channels. The 
solid lines show the bounds calculated directly through the 
MCMC analysis, while the shaded regions represent the 
ranges as obtained from Eq. ([1]) and Table [1] Indeed, the 
solid lines are fully covered by the shaded regions, as it 
should be. The vertical gray stripe indicates the range of 
WIMP masses (ttt-dm = 6 — 8 GeV) th at provide a g ood fit 
to CoGeNT and DAM A/LIBRA data (jHooper et al.|[20Tol) . 
The vertical dotted line marks the lowest DM particle mass 
5 GeV used in our PYTHIA simulations. This cut-off is not 
physical, but occurs because PYTHIA does not work be- 
low that energy. Therefore, the extrapolations of our results 
are shown below 5 GeV DM mass. Note that directly cal- 
culated lines for the fj,- and e-channel have slightly steeper 
slopes than given by the shaded regions, which increase as 
oc TTiDM- As the /-parameters generally fall off more slowly 
for lower ttt-dm, this behavior is also typical of the other 
channels. Thus, for lower ttidm values, one should actually 
get slightly stronger bounds than calculated directly from 
Eq. ([T}, and so our extrapolations shown in Fig. [S] are some- 
what conservative. 

We see that for WMAP 7, depending on the annihilation 
model, the limiting DM particle mass below which the up- 
per bound on the annihilation cross section drops below the 
standard thermal production value is in the range 4.5 — 10 
GeV, while the corresponding numbers reachable for the 
Planck and CVL experiment are 19 — 43 GeV and (45 — 100) 
GeV, respectively. Since the /.t-channel represents our most 
conservative case, one can instead say that, according to 
currently available CMB data, the annihilation cross sec- 
tion should be below the standard value of 3 x lO-^^ cm^s^^ 
as long as todm < 5 GeV. E3 Thus, for the CoGeNT and 
DAM A/LIBRA best-fit mass region 6-8 GeV, the stan- 
dard thermal production cross section is still compatible 



The numbers given here correspond to the l-sigma upper 
bounds. 
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Fig. 5. WMAP, Planck, and CVL l-sigma constraints on 
the (aAv) — Todm plane for /i (upper panel) and e (lower 
panel) annihilation channels. The solid lines show the up- 
per bounds on annihilation cross section as determined di- 
rectly through full MCMC calculations. The shaded re- 
gions around solid lines show the results from the sim- 
ple recipe of Eq. ([T]) with values of r taken from Table [TJ 
The vertical gray stripe shows the range of WIMP masses 
(rriDM = 6 — 8 GeV) tha t provide a good fi t to CoGeNT 
and DAM A/LIBRA data (iHooper et al.ll2010D . The vertical 
dotted line marks the lowest DM particle mass 5 GeV avail- 
able for PYTHIA simulations. Extrapolations are shown 
below that value. 



with the CMB measurements. This could quite possibly be 
changed soon as Planck results become available. Of course 
one should keep in mind that the cross section bounds given 
here directly apply to redshifts of 2 ~ 1000, and if the cross 
section depends on velocity (as in the case of P-wave anni- 
hilation or Sommerfeld-enhanced scenario), one should be 
careful in converting these numbers to the values relevant 
at 2 = 0. Also, note that the standard annihilation cross 
section 3 x lO-^^ cm'^s"^, with respect to what we are com- 
paring our CMB bounds to, provides the desired thermal 
relic density (i.e., rioM ~ 0.3) only if WIMPs are annihi- 
lating through S-wave processes. 

Except for a typical assumption of S-wave anni- 
hilation our results are largely mo del- independent, 
so all the particle-physics scenarios motivat ed by 
DAMA/LIBRA and CoGeN T results iFTtzpatrick et al 
2010; 
2010: 



Andreas et al 



Barger et al 



I2010al : IChne et al.l 



I2OIOI: IChang et all 120101: 
2010hl:" 



Hooper et al.l 



■Foot 
20101: 



Fitzpat rick fc Zurekl 120101 : lEssig et al. 1 120101: iBarger et all 



iKuflik et all I2OIOI: iFornengo et"all 



20 lit iBucklev et all 120111 )" including 

models that 
2OIOI: 



ahysics 



theoretically well-motivated particle pnysi 
predict light DM, s u ch as the MSSM (pjldman et al 



Bottino et al.l 
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201(| and the NMSSM (iKang et alJ 120 lit iBelikov et alJ 
20101: iGunion et all 120101: iDraper et al.ll201lD . are stri~ 



gently tested by WM AP and Plan c k. 

A recent study bv lDobler et al.l (|201lD suggests that an- 
nihilating DM with masses in the range jtidm = 1 — 20 
GeV (i.e. compatible with CoGeNT and DAMA best-fit 
region), along with annihilation cross section (ctaw) ~ 
9 X 10~^^ cm'^s"^ (i.e. 30 times the standard thermal relic 
cross section), may give a good match to the observed 
Fermi and WMAP haze. From Fig. [S] it is clear that, 
for the DM particle masses in this range, such high an- 
nihilation cross sections certainly conflict with available 
CMB data. If the required boost by a facto r of ^ 30 is 
obtain e d via the Sommerfeld enhancement ([Cirelli et al.l 
I2009bt lArkani-Hamed et al.l [20091: ISlatveij[2010[ ). then at 
z ~ 1000 the expected cross section is at least as large as 
the value given above, so the mod el could be in even big- 
ger tr ouble. The alternative cl aims |Goodenough fc Hoopeij 
120091: iHooper fc Goodenoughll201ll) of Ught DM annihila- 
tions in our Galaxy will also be stringently tested. 

Fina lly, it would be in teresting to compare the bounds 
given m ISlatver et all (^9) , which are based on the results 
obtai ned bvlGalli et al.. (|2009i) . with our values. However, 
smce iGalli et al.l (l2009l) used 'on the spot' approximation 
for the energy deposition, the comparison can only be 
approxima te. The WMAP5 2- sigma bound as given by 
Eq. (6) of ISlatver etaP (l2009l) can be cast in the form 
{pAv) [3 X 10-26 cm^s-i] /moM [GeV] < 0.12//. To com- 
pare this with our WMAP7 results, we have to set f,y — 
in Eq. ([T]) and, in the above relation, use a typical value 
for the /-parameter at z 1000, which from Fig. 0] is 
/ ^ 0.85. Thus, the value for r that should be directly 
comparable to the 2-sigma WMAP value given in Table [T] 
is 0.12/0.85 ~ 0.14. Considering the differences in the treat- 
ment for the energy deposition, this value agrees reasonably 
well with our result. However, there are significantly greater 
differences if one compares the forecasts for the Planck 
and CVL experiments. In our case Planck and CVL would 
tighten the 2-sigma bound of WMAP by a factor of '-^ 6 and 
~ 14, respectively. The c orresponding numb ers (~ 13 and 
~ 40, respectively) from iGalli et aL (j2009t ) are certainly 
more optimistic. For CVL, some of this d iscrepancy i s surel y 
due to the higher ^max value assumed in lGalli et al.l (j2009tl : 
^max — 2500 compared to our ^max = 2000. Even though 
we are comparing here the results derived from WMAP5 
and WMAP7, the difference is expected to be negligible 
because the accuracy of the cosmological parameters with 
an additional two years of WMAP data improves only very 
moderately (typically somewhere between 5 — 15%), and 
thus cannot be the cause of the above discrepancy. 

NOTE: After th e first version of th is work was submit- 
ted, a new paper bv lGalli et al.l (|201lt) appeared where the 
authors have performed a similar analysis (now also using 
WMAP7 data), this time finding results that agree with 
ours more closely: typical deviations are now within fac- 
tors 1.2 — 2, with their bounds on annihilation cross section 
being somewhat stronger. These relatively small deviations 
are actually quite negligible keeping in mind several ap- 
proximations used in both analyses. 

To see how much the omission of the additional excita- 
tions of atoms could change our results, we also performed 
several calculations where the exc iting Ly-g ph otons are 
treat ed along the lines presented in lChen fc Kam ionkowskil 
(|2004l') . In agreement with the claim in Galli et al.i ()2011f) . 



we find that the bounds on annihilation cross section are 
getting tighter only up to ~ 10%; i.e., the excitations seem 
to have only a relatively weak effect with the current pre- 
scription. 

Note on protons 

Although for most of the channels and, in particular, for 
the lower end of the considered DM particle mass range, 
the number of protons produced is completely negligible, for 
quark and gluon channels the contribution can reach ^ 15% 
of the total energy input if one has masses at the higher end 
of the considered range. As this is still only a relatively mild 
contribution, we did not attempt any detailed modeling 
for the proton component and simply assumed that all of 
this energy is absorbed as heat by the cosmic medium. A 
simple justification is the following. The energy loss rate 
for protons with the energies of interest in this paper is 
dominated by proton-proton scattering. Thus the loss rate 
r = —-^^ — npCCTppKpp, where Hp is the number density 
of target protons, app the scattering cross section, and Kpp 
the inelasticity parameter. The cross section tTpp depends 
only weakly on proton energy with typical values of 30 



40 mbarn and inelasticity parameter 



PPr 



0.5 for the 

energies of interest in this work (see e.g. IWommer et al.l 
2008). At redshifts z - 1000 we therefore get F ~ lO'^^ g-i 
for the energy loss rate. Comparing this to the expansion 
rate H at the same redshift, H ~ 0.4 x 10""'^^ s"""^, we see 
that F > iJ, so one might expect a significant fraction of 
the proton energy to be absorbed. 



5. Summary 

In this paper we have calculated the existing and future 
CMB constraints on annihilating DM assuming the S-wave 
annihilation mode to be dominant. Our results can be sum- 
marized as follows. 

• In full ag reement with our earlier findings, first pre- 
sented in iHiitsi et al.l (|2009t ). we confirm that in case 
of any realistic halo mass-concentration relation, along 
with annihilation cross sections {aAv} that do not vio- 
late the CMB data, the contribution from structure for- 
mation to the CMB signal is completely negligible. Also, 
the contribution to the low-redshift ionization fraction 
is very mild, and thus the annihilating DM models 
compatible with CMB measurements could only play 
a minor role in helping to reionize the low- 2 Universe. 
Therefore for low-mass DM the CMB constraints can 
be considered almost free of cosmological uncertainties. 

• At the same time, the particle physics uncertainties 
on the DM annihilation channels can all be described 
with one single parameter /^, the energy fraction car- 
ried away by neutrinos. Therefore the DM annihilations 
to two muons represents the least stringently tested DM 
scenario. 

• As our main results, in the form of Eq. ([T]) and Table [U 
we provided a simple recipe for estimating the upper 
bounds on annihilation cross section, valid for a broad 
range of DM models. Two examples for 1-sigma upper 
bounds for the fi- and e-channels (which are the two ex- 
tremal cases among all of the channels considered in this 
paper; i.e., those upper bounds bracket the expected re- 
sults for more general annihilating DM models, given 
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as a superposition over the basis channels) are shown in 
Fig.H 

• For the DM particle masses todm = 6 — 8 GeV, which 
give best fits to CoGeNT and DAMA/LIBRA data, cur- 
rent CMB data is still compatible with a standard ther- 
mal relic annihilation cross section 3 x 10^^^ cm'^s^^ 
only if the annihilations are dominantly into the /i or 
r (and corresponding 4 lepton) channels. All other an- 
nihilation channels already now conflict with the CMB 
data. 

• The sensitivity of Planck space mission allows one to 
test all the DM annihilation channels definitively. If the 
Planck mission will not find signals of annihilating DM 
in the CMB spectrum, all the light DM scenarios mo- 
tivated by DAMA/LIBRA, CoGeNT, Fermi haze, and 
WMAP haze will be stringently tested if DM is a ther- 
mal relic. 

• Our findings strongly disfavor a claim that annihilating 
DM with moM 1 - 20 GeV and {(jav) -Ox 10-^5 
cm'^s"^ (i.e. ~ 30 times above the standard value) could 
be a cause for Fermi and WMAP haze. 

• The last conclusion applies to all DM scenarios with 
large boost factors from the Sommerfeld enhancement. 
Because CMB is sensitive to DM annihilations at z ^ 
1000, at that time the Sommerfeld enhancement of anni- 
hilation cross section should have been at least as large 
as today. Therefore, for light DM scenarios the boost of 
DM annihilation in our Galaxy can only come form the 
halo substructure. 

• Our MCMC calculations assumed standard 6-parameter 
ACDM cosmology, which was extended by two addi- 
tional parameters describing the annihilating DM: mDM 
and {aAv) I'm-DM- It turns out that, after introducing 
these 2 additional degrees of freedom, all the previous 
6 parameters, albeit with small shifts (the most notice- 
able of which being the ones for Aa, and ns), were still 
as precisely determined. 

• The weakest point of the analysis presented in this pa- 
per is the frequently used simple approximation about 
how the input energy of the cascading particles gets 
partitioned between ionizations and heating of the envi- 
ronment. For better treatment of this issue one possibly 
has to rely on detailed Monte Carlo calculations. This, 
however, we leave for a possible future study. 
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